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LINEAR CHAOS FOR THE QUICK-THINKING-DRIVER 

MODEL 
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Abstract. In recent years, the topic of car-following has experimented 
an increased importance in traffic engineering and safety research. This 
has become a very interesting topic because of the development of driver¬ 
less cars [23] . 

Driving models which describe the interaction between adjacent ve¬ 
hicles in the same lane have a big interest in simulation modeling, such 
as the Quick-Thinking-Driver model. A non-linear version of it can be 
given using the logistic map, and then chaos appears. We show that an 
infinite-dimensional version of the linear model presents a chaotic be¬ 
haviour using the same approach as for studying chaos of death models 
of cell growth. 


1. Introduction 

Dynamical systems on infinite-dimensional linear spaces allow us to de¬ 
scribe phenomena in which a chaotic behaviour is exhibited. The dynam¬ 
ics of the solution Co-semigroups of gene amplification-deamplification pro¬ 
cesses with cell proliferation have been widely studied by Banasiak et al. 
[21 0 III E] (see, also, [D m EZ1). Such a behaviour can also be found on 
certain size structured cell populations, c.f. [201 EH- 

These models are very similar to the ones used to describe traffic. Models 
of car-following describe the interaction between adjacent vehicles in the 
same lane and show how one vehicle follows one another. Since the studies of 
Greenshields [25, 26] in the 1930’s, several models started to appear sharing 
this point. See for instance the work of Pipes [33], Helly [29], and Herman, 
Montroll et al mm in the 1950’s. These models were perefectioned and 
improved. For a historical evolution of these models we refer the interested 
reader to m- 

During the 2000’s, great advances in robotics and artificial intelligence 
have permited to develop prototypes of autonomous cars. The advances can 
be followed looking at the DARPA Grand Challenge competition [16] , In 
fact, motor companies like Audi or BMW have presented their first models 
very recently m- 
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We would like to emphasize that these seminal models, which are based 
on simple linear equations, cannot describe all the highly complex situations 
that, actually, occur on a roadway. These theoretical models work better 
on long stretches of road with dense traffic. In this note we consider one of 
these models, the Quick-Thinking-Driver model, for an infinite number of 
cars on a track. The infinite system of ordinary differential equations that 
are required to model the behaviour of these vehicles can be represented as 
a linear operator on an infinite-dimensional separable Banach space. Then, 
using some results of linear dynamics of Co-semigroups we can prove the 
existence of different chaotic behaviours for the solutions of these equations. 

2. Preliminaries 

In the basic formulation of any of the car-following models, there is a 
relation between the acceleration of a car and the difference between its 
velocity and the velocity of the car that goes in front of it. We assume that 
a driver adjusts her speed according to her relative speed respect to that 
car. This can be described by means of the following equation 

(1) u[(t+ ti) = \i(u 2 {t) - 

in which u\ stands for the velocity of the car 1, U 2 for the velocity of the 
car in front of the car 1, namely car 2 , t\ denotes the reaction time of driver 
1 , and the positive number Ai is a sensitivity coefficient that measures how 
strong the driver 1 responds to the acceleration of the car in front of her. 
Usually Ai lies between 0.3 — 0.4s' - 1 m- Under the assumption that all 
drivers react “very quickly”, one can take t\ = 0. This is known as the 
Quick-Thinking-Driver (QTD) model. 

If we consider this model for an infinite number of cars circulating on a 
road, then the QTD model is given by an infinite system of first-order ordi¬ 
nary differential equations. As far as we know, this has not been considered 
yet in the literature: 

Definition 2.1 (The Infinite Quick-Thinking-Driver model). Let ( Ui(t))i 
denote the velocities of an infinite number of cars on a track. Suppose that 
for every i 6 ff, the car i behaves following the (QTD) model respect to 
one in front of it, the car i + 1. The following infinite system of ordinary 
differential equations represents this situation. 

(2) u'fit) = \i(u i+ i(t ) - Ui(t)) for i <E N. 

We also assume that the velocities at t = 0, («j(0))j, are given. We call 
to this model the Infinite Quick-Thinking-Driver model. We will refer to it 
as the (IQTD) model. 

McCartney and Gibson studied chaos for the (QTD) model making the 
parameter Aj to be proportional to the velocity (( 2 j) , i.e. \ is replaced by 
7 Ui(t), with 7 being the new sensitivity coefficient, and making the velocity 
of the leading car to be constant. This converts dTJ) into the logistic equation, 
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see also [301131] and then chaotic phenomena can be found in its solutions. 
We will see that even in the linear (IQTD) model given by ([2]) chaos still 
appears. 

Let X be a separable infinite-dimensional Banach space. We assume that 
the reader is familiar with the terminology of Co-semigroups on Banach 
spaces, see for instance [22 ]. 

Next, we recall some basic definitions on linear dynamics of Co-semigroups. 
A Co-senrigroup {T t }t> o on X is said to be hypercyclic if there exists x € X 
such that the set {T)x : t > 0} is dense in X. An element x € X is a 
periodic point for the semigroup if there exists t > 0 such that Ttx = x. A 
semigroup {T t } t > o is called Devaney chaotic if it is hypercyclic and the set 
of periodic points is dense in A. In this setting, Devaney chaos yields the 
sensitive dependence on the initial conditions, as it was seen by Banks et al 

mm- 

Another extended definition of chaos is distributional chaos [32]. A Co¬ 
semigroup {Tt}t> o on X is said to be distributionally chaotic if there exists 
an uncountable subset S C X and 5 > 0 such that, for each pair of distinct 
points x,y G S and for every e > 0, we have Dens({s > 0; \\T s x — T s y || > 
5}) = 1 and Dens({s > 0; ||T s x — T s y || < e}) = 1, where the upper density 
of a set A C R is defined as 

fi(An [0,t]) 


Dens(A) = limsup 

t —>oo l* 

where n denotes the Lebesgue measure on M. We will show that both De¬ 
vaney and distributional chaos appear in the solutions of the (IQTD) model 
using a suitable underlying Banach space. 

3. Chaos for the Quick-Thinking-Driver model 

Let us consider C(s), 0 < s < 1, the weighted space of summable se¬ 
quences defined as 


(3) 4(s) = < (ui)teN € : IK^ieNlls = ^ < oo > 

l «eN J 

If s < 1, then we can consider the velocities of all the cars of the (IQTD) 
model in i\ (s). We point out that such a choice of weights gives more 
importance to velocities of cars with low index i. 

In order to solve the system of equations in ([2]), we pose the following 
abstract Cauchy problem on ^i(s) 


(4) 


u'{t) = Au(t ) 

u( 0) = (u;(0)); e N 


where the operator A is defined as 

(5) (Au(t))i = Xi(u i+ i(t) - Ui(t)), i € N 
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for u = G 4(s) and t > 0, being (uj(0))igN the vector of velocities 

of the cars at t = 0. 

If the supj eN | A,; | < oo, then the solution to (j4]) can be represented by a Co¬ 
semigroup {Tt]t> o on ^ ,1 ( s ) whose infinitesimal generator is A with domain 
the whole space D(s) . Then the operators in this Co-semigroup can be 
represented as T t = e tA = Y1T= o(^4) n / n - for all t > 0, see for instance 
[22]*Ch. I, Prop. 3.5. 

When considering in the death models used for modeling cell-growth 

are of the form 

(6) u'fit ) = j3iUi+i{t) — a.iUi(t ) with >0 for i G N 

In proposition 7.36 of m they proved the following result: 

Theorem 3.1. If sup,- a t < liminfj^oo fy, then the solution Co-semigroup 
of (|6|) is Devaney chaotic (and topologically mixing) on I i(l). 

More details can be found in mm- This proof is based on an application 
of the Desch-Schappacher-Webb (DSW) criterion [18], see also [5; ‘IQj 122] for 
other reformulations of this result. The application of the (DSW) criterion 
relies on verifying that the point spectrum of the infinitesimal generator of 
the Co-semigroup contains “enough” eigenvalues. A criterion directly stated 
in these terms was firstly given for operators by Godefroy and Shapiro in 
[SI]• We will use the following version of the (DSW) Criterion, see [2?j*Th. 
7.30. 

Theorem 3.2. Let X be a complex separable Banach space, and {T t }t> o 
a Co-semigroup on X with infinitesimal generator (A, D(A)). Assume that 
there exists an open connected subset U and a weakly holomorphic function 
f : U —> X, such that 

(1) UniR ± 0, 

(2) /(A) G ker(A I — A) for every A G U , 

(3) for any x* G X* , if (/( X),x*) = 0 for all A G U, then x* = 0. 

Then the semigroup {T t }t> o is Devaney chaotic (and topologically mixing). 

Since we can assume 0.3 < Aj < 0.4 for all i € N, c.f. m, we can apply 
this criterion to the (IQTD) model. The proof of the following result can be 
compared with the one of [27]*Prop. 7.36. 

Theorem 3.3. Let 0 < a < (5. Let us also consider the (IQTD) model on 
with 0 < s < a/fi. If a < A* < /3 for all i € N, then the solution 
Co-semigroup of the (IQTD) model is Devaney chaotic on £\(s). 

Proof. Fix s with 0 < s < a//3 and consider the solution Co-semigroup of 
the (IQTD) model whose infinitesimal generator is given by the operator A 
defined in ([5]). For applying Theorem 13.21 we consider as U the open disk 
centered at 0 and radius 0 < e such that (e + (3)s/a < 1, i.e. U = {p, G 
C : \/a\ < e}, and then (1) of Theorem 13.21 holds. We define the function 
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/ : U —> £i(s) as /(/j) = h^, with h^ defined as 



Then f(/j) £ £i(s) because, by the choice of we have 
( 8 ) 




i—1 


(9) 


< oo. 


In fact, f(/u) is an eigenvector of A associated to the eigenvalue /i. This 
gives (2) of Theorem 13.21 Finally, for proving (3) of Theorem 13.21 let 7 £ 
^ 00 (1/s)j the dual space of Q(s), such that (/(^), 7 ) = 0 for all A £ U. The 
function p —> (/(^), 7 ) defined as 



( 10 ) 


converges uniformly for all /a £ U by[U which gives that / is weakly holo- 
morphic on U. 

If we consider p, = — X± £ U we have that (flOl) is reduce to 71 , which gives 
7 ! = 0. Then, we divide the expression in (fTOl) by p + \i- The result is still a 
weakly holomorphic function on U. Now, if we consider p = —A 2 £ U, then 
we get 72 = 0. Repeating the process inductively, we get that all 7 j = 0, 
i £ N, and the proof is finished. 


□ 


Clearly, we get Devaney chaos for the aforementioned experimental values 
obtained for the (IQYD) model. 

Corollary 3.4. Let us consider the (IQTD) model on ^i(s), with 0 < s < 
3/4. If 0.3 < A i < 0.4 for all i £ N, then the solution Co-semigroup of the 
(IQTD) model is Devaney chaotic on t\ (s). 

As a result, if an infinite numbers of cars drive along a road and all 
the drivers respond in the same way to any change in the velocity of the 
car in front, this may lead to a chaotic situation, depending on the initial 
velocities of the cars. Therefore other considerations should be taken into 
account in order to manage the speed adjustment. This also contrasts with 
the behaviour of the Co-semigroup solution on t\ studied in m where it is 
proved that this model is not Devaney chaotic on C(l). 

Up to here, we have studied Devaney chaos for the (IQTD) model. It is 
also well known that distributional chaos holds whenever the DSW criterion 
can be applied BW As a result, we obtain the following corollary. 
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Corollary 3.5. Let us consider the (IQTD) model on ^i(s), with 0 < s < 
3/4. If 0.3 < A i < 0.4 for all i £ N, then the solution Co-semigroup of the 
(IQTD) model is distributionally chaotic on i\ (s). 

Roughly speaking, this means that we can pick two vectors of initial ve¬ 
locities from an uncountable set such that there will be intervals of arbitrary 
time length in which the velocities of the cars are very similar and intervals in 
which there exists at least some positive difference between their velocities. 

We conclude this section with a comment regarding the stability of the 
solutions. An analysis of local and asymptotic stability of the (QTD) model 
was considered in [28) . From the point of view of Co-semigroups, we recall 
that a Co-semigroup of the form {e tA }t> o defined on a Banach space X is 
(uniformly) exponentially stable , [ 22 ], if there exists rj > 0 such that 

( 11 ) lim e^Ne^ll = 0 . 

£—>■00 

This situation is not fulfilled in our case. Nevertheless, a weaker version of 
this condition can also be considered. We say that {e tA } t >o is exponentially 
stable on a subspace Y C A if there exists > 0 such that for any y £ Y 
we have 

( 12 ) lim e^||e tA y|| = 0 , 

£—>•00 

This analysis is sometimes considered when studying chaos of Co-semigroups, 
as it is done in pouees]- 

Theorem 3.6. Let 0 < a < f3. Let us consider the (IQTD) model on 
£i(s), with 0 < s < a/(5. If a < A* < @ for all i £ N, then the solution 
Co-semigroup of the (IQTD) model is exponentially stable on the subspaces 

(13) Y s := span {h^ : \fj,\ < e, 5ft(/z) < <5}, 
for every — e < 5 < 0, where e > 0 satisfies s(e + (3)/a < 1. 

Proof. Fix 0 < s < a/(5, £ > 0 satisfying s(e + (3)/a < 1 and — e < 5 < 0. 
Take 0 < p < —5 and y £ Ys of the form y = Yli =i a il l in ■ Define 6 y = 
max{3?(/Xj) : 1 < i < k}. Clearly, r] + d y < 0, then 



which tends to 0 when t tends to oo. □ 

Remark 3.7. For the experimental values of 0.3 < A * < 0.4 for all i € N, 
Theorem l3.6l states that there are subspaces where the solution Co-semigroup 
of the (IQTD) model is exponentially stable on C(s), for 0 < s < 3/4. 
Moreover, in these spaces there are also subspaces where the solution Co- 
semigroup of the (IQTD) model is not exponentially stable. If e > 0 satisfies 
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s(e + (3)/a < 1, we can take, for instance, the case where initial conditions 
are taken on Z := spanj/i^ : \/j,\ < e, > 0}, since we have that 

lim^oo HT^IIs = oo for all z € Z \ {0}. 

4. Final comments 

In this note we have related the dynamics of models of traffic following 
with the dynamics of death process of cell growth. We also point out that 
some other models like the one of forward and backward control j28j can 
also be related with birth-and-death process of cell-proliferation. 

It is fair to note that the chaotic properties reported in this paper refer 
to the whole-space of solutions whereas only positive solutions make sense. 
Therefore, being important, e.g., in the analysis of the stability of numeri¬ 
cal schemes, does not necessarily say anything important about the actual 
behaviour of the “physical” trajectories of the cars, in the same way as it 
is commented in [4j. To sum up, in dense traffic drivers follow one another 
very closely and small disturbances such as acceleration or deceleration of 
one vehicle might be passed over or amplified along the line of vehicles on 
the road. These disturbances can be considered in ((4]) by adding to the right 
part of the equations a term related to the notion of acceleration noise (a 
kind of measure for the smoothness or the quality of traffic flow, see m- 
A chaotic flow appears and this can lead to an impredictible dynamical 
behaviour and, in some cases, to accidents. 
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